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Motivated by the possible charge ordering in the recently discovered superconductor 
Na2,Co02 yH20, which at filling x — 1/3 would correspond to a half-filled honeycomb lattice, 
we investigate the single-hole dynamics of the t-J model on this lattice. Using exact diagonaliza- 
tion, series expansion and the self-consistent Born approximation, we calculate the quasi-particle 
dispersion, bandwidth and residues and compare our findings with the well established results for 
the square lattice. Given the similarities between both lattices, it is not surprising that we find 
qualitatively similar features for the honeycomb lattice, namely an almost flat band along the edge 
of the magnetic Brillouin zone and well deflned hole pockets around the corners. However, we report 
a considerable disagreement between the three methods used concerning the bandwidth and discuss 
possible origins of this discrepancy. 
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I. INTRODUCTION 

The study of doped antiferromagnets has been a cen- 
tral subject in the area of strongly correlated electron 
systems. The pioneering work by Bulaevskii, Nagaev and 
Khomskiii as well as Brinkman and RiceS, led the way to 
innumerable works, especially motivated by the discovery 
of high-temperature superconductors'^ and the suggestion 
by Anderson^ that the physics of these materials could 
be described by a two-dimensional Hubbard model on the 
square lattice with large on-site Coulomb repulsion. In 
this limit, the Hubbard model can be mapped to the t-J 
model^, which has been shown to describe experimentally 
determined single-hole properties very accurately, consid- 
ering t, t' , t" and J as independent variablea^iSiS. At 
half-filling, the t-J model reduces to a Heisenberg anti- 
ferromagnet, because of the single occupancy restriction. 
The ground state of the Heisenberg model on a bipar- 
tite lattice is the quantum Neel state with long-range 
magnetic order. Doping away from half-filling allows the 
charge carriers to propagate and hence scramble the Neel 
order, leaving strings of overturned spins behind. A va- 
riety of tools has been employed to study the hole mo- 
tion in an antiferromagnetic backgrounds^, e.g., exact di- 
agonalizations (ED)^^i^^i^'^i^"^i^^i^^i^''', the self-consistent 
Born approximation (SCBA)i2ii2i22i2i, quantum Monte- 
Carlo simulation32242i24 and series expansion^^. In the 
parameter range relevant for the copper-oxide materials, 
i.e., J/t w 0.3, these methods agree reasonably well2^ 
and the main features emerging are a minimum of the 
quasi-particle dispersion at (§,§), with a very flat band 
along the edges of the Brillouin zone and a practically 
vanishing quasi-particle residue at (tt, tt). 

The recent discovery of superconductivity in 
Na£cCo02 yH20^® has heated interest in models on 
the triangular lattice, to which the SCBA has already 
been applied22i2£. In this material, the two-dimensional 
C0O2 planes are separated by insulating layers of 
Na ions and water molecules. The cobalt ions at 



the center of the CO2 octahedra form a triangular 
lattice. Superconductivity has been observed for 
doping 0.22 < x < OAT^^i^^. Spin and charge 
ordering tendencies in Naa;Co02 have been addressed 
recentl y?^'??i?^i?^ , and could, at flUing x = 1/3, lead 
to an effectively half-flUed honeycomb lattice, where 
hole and electron doping is possible. Estimations for 
the parameters in the t-J model converge at around 
J/t sa 0.1 — 0.32s., slightly smaller than for the copper 
oxides. 

Motivated by the possible application to this new su- 
perconductor, we investigate the dynamics of a single 
hole on the honeycomb lattice, using exact diagonaliza- 
tion, series expansion and the SCBA and compare our 
findings with well established results for the square lat- 
tice. The rest of the paper is organized as follows. In 
Sec. m we define the model and introduce the notations 
for the lattice geometries under consideration. Sec. IIIII 
then contains a summary of the most important features 
of the applied methods. We first consider the simplified 
t-J^ model in Sec. II VI with Ising rather than Heisenberg 
spin interactions where quantum fluctuations are absent 
and present our main results for the full t-J model in 
Sec. 13 Finally, Sec. IVII contains our conclusions. 



II. MODEL 



The t-J model originates from an effective low-energy 
description of the Hubbard model in the limit of strong 
Coulomb repulsion {U 3> t), where double occupancy is 
forbidden and J = At^ /U . An accurate treatment of the 
large-f/ limit of the Hubbard model^ leads to correlated 
hopping terms, which we neglect in this work, in order to 



simplify the discussion. The Hamiltonian therefore reads 



<i,J>, cr 



(1) 



Here c^^ are the usual electron creation operators, with 
Tii — c|^|Ci| + cl^Cii, Si are the electron spin operators 
and P the projection operator that excludes double oc- 
cupancy. The square (honeycomb) lattice has coordina- 
tion number Z = A (Z = 3). To simplify comparison 
between exact diagonalization, series expansion and the 
SCBA, we assign the zero energy level to the state with 
one static hole 



Eo - i;H"-">^-g + J/2 = 



(2) 



Here C takes into account the quantum fluctuations. Ac- 
cording to series expansion, ^ — 1.08 for the square lat- 
ticeSSi, and C = 1-10 for the honeycomb lattice, whereas 
according to exact diagonalizations on the honeycomb 
lattice, we find C ~ 1-182 for TV = 18 and 1.165 for 
N = 24. 

In order to gain additional understanding of the accu- 
racy of the applied methods and the mechanism of hole 
propagation, it is insightful to study the simplified t-J^ 
model 



<i,j>, a 



<ij> ^ 

with Ising rather than Heisenberg spin interactions. In 
the absence of spin fluctuations, C = 1 in Eq. 0. 

We study the above models on the square and the 
honeycomb lattice. The latter one being not a Bra- 
vais lattice, is best described as a triangular lattice with 
two sites per unit cell. The reciprocal lattice is then 
again triangular, and the first Brillouin zone is a honey- 
comb, as shown in Fig. ^ Highly symmetric points are 
r= (0,0), X= (f ,0), L= (f ,0), y = ('^'^)' and 

W = ^ i^^) ' Si^ila'^ drawings for the square lattice 
are shown in Fig. [21 with the main purpose to intro- 
duce the notation for points of special interest F = (0, 0), 
M = (7r,0), X = (7r,7r), and S ^ (f , f ). 

III. METHODS 

We calculate several single-hole properties, namely 
the quasi-particle dispersion, the bandwidth and the 
residues, using exact diagonalization and series expan- 
sion methods as well as the SCBA. We briefly recall the 
ideas, advantages and possible difficulties associated with 
these three methods. 
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FIG. 1: Honeycomb lattice (described as a triangular lattice 
with two sites per unit cell) in real (a) and reciprocal space 
(b). The magnetic (shaded) and the full Brillouin zone coin- 
cide. The labels are defined in the text. Note that in contrast 
to the square lattice, the white area does not correspond to 
the Brillouin zone, but to the periodicity of the quasi-particle 
weight measured in ARPES experiments. 
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(b) 



FIG. 2: Square lattice in real (a) and reciprocal space (b). 
The shaded area corresponds to the magnetic Brillouin zone. 
Highly symmetric points are mentioned in the text. 



A. Exact diagonalization 

In exact diagonalization one can rather easily calcu- 
late spectral functions directly in real frequency, using 
the so called continued-fraction techniqueiS*^. This al- 
lows to obtain unbiased results, including residues, al- 
though only on finite-size samples. In special cases it 
is possible to calculate single-hole spectral functions for 
t-J models on samples up to 32 sites, depending on the 
lattice symmetries. This method has been applied on nu- 
merous occasions to obtain single hole spectral functions, 
especially for the square latticeii^iSiiiMiiSiiSiii. To our 
knowledge single-hole properties for the honeycomb lat- 
tice have not been studied in ED before. In the following 
we will mainly present results obtained on iV = 18 and 
= 24 sites samples. They share the advantage of hav- 
ing the symmetries of the infinite lattice and furthermore 
have the W point in their discrete Brillouin zone, there- 
fore giving a reasonable estimate of the bandwidth. For 
some quantities we considered systems up to A^ = 32 
with one hole. 
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B. Series expansions 

The Ising series expansion for the single-hole disper- 
sion of the t-J model has been done previously for the 
square latticeS^. To perform the series expansion, one 
needs to introduce an expansion parameter x, and divide 
the Hamiltonian into an unperturbed part (i?o) ^^^id a 
perturbation {V) as follows 

H/J = Ho + xV , 

with 

(ij) i 

with a staggered z-field of strength r to improve conver- 
gence. This term vanishes in the total Hamiltonian in 
the isotropic exchange limit x —^ 1. The Ising part of the 
exchange, Hq, forms the unperturbed Hamiltonian while 
the perturbation V includes both the hole hopping terms 
and transverse spin fluctuations. The original Hamilto- 
nian for the t-J model is recovered in the limit x = 1. We 
expect the series expansion to give more reliable results 
in the small t/ J region, where the perturbation V is less 
important. The series for the single- hole dispersion has 
been computed to order x^^ for both the square and the 
honeycomb lattice, involving a list of 325737 and 28811 
clusters with up to 13 sites, respectively. This extends 
the previous series for square lattic^^ by only 2 terms. 
We also compute the new series for the quasi-particle 
residue, again to order x^^. 

For the t-J^ model, one can also form a direct expan- 
sion in powers of t/ J^. This series has been computed 
to order {t/J^Y^ [(t/J^)^^] for square [honeycomb] lat- 
tice requiring a hst of 4654284 [80905] clusters with up 
to 15 [14] sites. Previous results for the t-J'^ model on 
the square lattice^^ used expansions of order (t/J^)^^. 

These series are available from the authors upon re- 
quest. 

C. Self-consistent Born approximation 

From a general point of view, the SCBA works well for 
a collinear state, whenever vertex corrections are small. 
In the case of the hole-spin-wave vertex, single loop cor- 
rections are zero, and higher order corrections are small 
numericall}i2i. Despite its simplicity, the SCBA has the 
disadvantage that its predictions are not immediately 
comparable to other methods or experimental results, be- 
cause it distinguishes pseudo-spin, i.e., holes living on a 
particular sublattice, rather than spin variables. 



1. Introduction 

Let us start this reminder by defining the bare hole 
operators d^, so that d\ cx c^i on the "f-sublattice {A) 
and d\ cx Cii on the |-sublattice {B). In momentum 
representation 

V jes 

N is number of sites and m = |(0|S'f|0)| the average sub- 
lattice magnetization, with to « 0.3 for the square lattice 
and TO ~ 0.24 for the honeycomb lattice. We denote the 
ground state of the undoped antiferromagnet, i.e. the 
quantum Neel state by |0). In the Ising case, there is no 
reduction of the sublattice magnetization due to quan- 
tum fluctuations and therefore to = 0.5. In this case, |0) 
is the classical Neel state. The quasi-momentum k is lim- 
ited to the magnetic Brillouin zone and the coefficients 
in ^ provide the correct normalization 

The retarded hole Green's function is defined as 

poo 

Gd.A{^,k) - -t / (0|dk,A(r)4^(0)|0)e"^dT , (5) 
Jo 

with the corresponding spectral function 

Ad.A(e,k) = Im Gd,yi(e,k) , 

TT 

and similar for Gd.si^i k). 

For spin excitations, the usual linear spin-wave theory 
is used^^ . It is convenient to have two types of spin waves, 
with — —1, and with — +1, and q restricted 
to the magnetic Brillouin zone. The spin operators can 
then be expressed as 

q 
q 

Although spin-wave theory is an expansion in powers of 
1/S, its application to the t-J model is restricted to the 
case S — 1/2, because only in this limit is it possible to 
flip a spin by absorption or emission of a single magnon. 
We therefore set 5 = ^ from now on. 
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Since the honeycomb lattice lacks inversion symmetry, 
7q and the parameters of the Bogoliubov transformation 
are complex: 

P 




where p is the vector to nearest neighbors and 5q is the 
phase of 7q. In this notation, the spin- wave dispersion is 
equal to eq = JZ jluj^. Hopping to a nearest neighbor in 
the Hamiltonian J^l gives an interaction of the hole with 
spin waves 



with 



iih,sw = ^ (5k.q4+q,.l^k,i3aq + 5k,q4+q,sCJk,A/3q + H.C^ , (6) 



k,q 



5k,q = (0|aqrfk,B|i?tMk+q,A|0) 



2t 



N(l/2 + m) 



E e-'-'-^+^C^+'i)- ((0|aqc]^c,T4^QT|0) + (Ojaq^ c.^cl^ c,t |0)) 



I E e-'^-'^+^C'+'i)- (0|aq {S- + S-) |0) ^ Zt^i^^u*^ + j^^^v^) , 



using the usual mean-field ground state factorization ap- 
proximation 

(0|aqc]^CjT4^c,t|0) « (0|aqcJ^c,T|0)(0|4^c,t|0) 
= (0|aqS'7|0)(l/2 + m) . 

For a coUinear state, the interaction © forbids single 
loop corrections to the hole-spin-wave vertex and it has 
been shown that two and higher order loop corrections, 
which correspond to Trugman processes^S are small2S*2i. 
This justifies the SCBA for Jmin < J/t <^ 1 according 
to which the hole Green's function satisfies the Dyson 
equation 

Gd{e, k) = - E l.9k-q.ql' Gd{e - Cq, k - q) -I- lO^ . 

(7) 

Clearly, there is a Jmin below which the ground state 
of the system is found in the sector \S^^^\ > 1, because 
in the limit J/t the Nagaoka theorem applies, pre- 
dicting a completely polarized ground state. A DMRG 
calculation^- for the square lattice with up to 9 x 9 sites 
estimates Jmin/t ~ 0.02 — 0.03. 

Due to the definition of the operators the Green's 



function jSJ is invariant under translation with the in- 
verse vector of the magnetic sublattice Q. 

G<i(e,k + Q) = G<i(e,k) . 

For the square lattice Q — (±7r, ±7r), and for the honey- 
comb lattice Q = (±^,±-2|). 

2. SCBA and ARPES experiments 

As explained in detail in Ref. |^ the imaginary part of 
the hole Green's function Gd is not directly measurable in 
ARPES experiments, because the incoming photon does 
not distinguish the sublattices. For comparison with ex- 
periments, we define an annihilation operator with given 
spin cr 

i 

The summation extends over all TV sites and k is the 
momentum of the free photoelectron detected in the ex- 
periment, i.e., it is not restricted to the Brillouin zone. 
The normalization is chosen to fulfill (0|c|^g.Ck(T|0) = 



5 



(a) 



k ^ k-q 

o 

(b) 



q 

^ k-q 

(c) 



FIG. 3: Vertices representing a single hole creation (a), a 
hole plus spin-wave creation (b) and the usual hole-spin-wave 
interaction (c). Bold lines correspond to Gc, solid lines rep- 
resent Gd, and dashed lines are spin waves. 
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FIG. 4: Dyson equation relating the hole Green's functions Gc 
(bold lines) and Gd (solid lines). Dashed lines represent spin 
waves. 



Si cLcicrlO) = 1. In this high-energy approxima- 
tion, every site is regarded as the possible origin of the 
emission of a free photoelectron. For the square lattice, 
this definition coincides with the normal Fourier trans- 
formation, but for the honeycomb lattice, this is not the 
case, since the full lattice is not a Bravais lattice. Eq. © 



J 



can therefore not be inverted in general. The spectral 
function measured in ARPES corresponds to the imagi- 
nary part of the retarded Green's function 



Gc(e,k) 



(0|cL(r)ck,.(0)|0)e'-dr. (9) 



The operator Ck.x acting on the quantum Neel state can 
either produce a single hole state, by removing an elec- 
tron with spin a —I from sublattice B, or it can produce 
a hole plus spin-wave state, by removing such an elec- 
tron from sublattice A. We adopt the notation of Ref. 
and denote the corresponding amplitudes by at and &k,q, 
respectively 



Ok = (0Mk,i3Ck,i|0) = V1/2 + ™ 

^k,q = (0|/3qdk-q,ACk,i|0) 

N{l/2 



These vertices are shown in Fig. |31 and afterwards are 
used in Fig. 0] to derive the Dyson equation for the 
Green's function of a hole with fixed spin Gc. In ana- 
lytical form, the Green's function Gc ® is given by 



Gc{e, k) = alGdie, k) + ^ |6k,q|' Grf(e - e^, k - q) + 2akGd(e, k) 



+Gd(e,k) 



^k,q5k-q,qGd(e " Cq, k - q) 



Re (&k,q.gk-q,q) Gd(e - Eq, k - q) 
Y ^k,q5k-q,qGd(e " Eq, k - q) 



L q 



(10) 



The numerical solution of ((7|l is straightforward. We re- 
place iO by 0.05* and use 32 x 32 points in the Brillouin 
zone and an energy interval Ae < 0.002 to carry out 
the iterations. Convergence is reached after around 50 
iterations for small values of J/t. Gc is subsequently cal- 
culated according to Eq. (^U)). Figs. 13 a) and (b) show 
the spectral functions Ad and Ac for a hole on the square 
lattice with J/t = 0.3, and on the honeycomb lattice 
with J/t— 0.2, respectively. A hole with fixed spin Gc 
has the same dispersion as a hole linked to a particular 
sublattice Gd, but clearly, the areas under the peaks is 
different. A discussion of the origin of these peaks and a 
comparison with ED calculations is presented in Sec. lVBl 
At this point, it seems useful to comment on the phys- 
ical difference between the two Green's functions. The 
electron creation operator Cj^ ^ does not correspond to a 
true quasi-particle of the system with long-range antifer- 
romagnetic order, because the spin a is not conserved. 
However, the hole creation operator associated with a 
particular sublattice dj^ ^ ,^ is well defined in this case 



and the quasi-particle residue corresponds to the area 
under the peak of the spectral function. In contrast, ex- 
act diagonalization and series expansion do not rely on 
the existence of long-range order and give direct access 
to the Green's function Gc. While there is no ambigu- 
ity concerning the single-hole dispersion, one has to be 
careful comparing residues, which we define as 

Z^{k}^ r A^(e,k)de, 

J —oo 

where Cp is the cut-off chosen to include only the lowest- 
lying peak. This definition has the disadvantage that 
is not well defined and that the tails of higher excita- 
tions are also taken into account. To circumvent these 
problems, we make the assumption of Lorentzian peaks 
and easily deduce the residue. Let us finally mention the 
finite-size effects that come into play in the calculation 
of Gc. While for the hole Green's function with fixed 
pseudo-spin, Gd, grids with more than 32 x 32 points in 
the Brillouin zone yield essentially the same results and 
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(b) Honeycomb lattice 



- A/E, k=r) 

- A (E, k=r) 




FIG. 5: (a) Spectral functions of a single hole in the t-J model 
on the square lattice with J/t — 0.3 at the minimum S and 
the maximum F of the dispersion, (b) Similar results for the 
honeycomb lattice with J/t = 0.2 at the minimum W and the 
maximum T of the dispersion. 



the energy grid with Ae = 0.002 fulfiUs the sum rule 
with good accuracy, i.e. the error is smaller than 0.5%. 
As can be deduced from the sum rulesSi, the situation is 
more subtle in the case of Gc 

/oo 
Ad (e, k) & = 1 , 
-OO 

J-oo q 

Clearly, a grid of only 32 x 32 points in the Brillouin 
zone does not lead to m = i — |uq|^ w 0.24, and 

the error in the sum rule is of around 3%. We solve 
this problem by choosing m according to the grid, i.e., 
~ 0.26. 

3. J/t limit 

Below Jmin/t, the ground state is found in the sector 
15^^(1 > 1, which cannot be accessed within the SCBA. In 
other words, in this limit, the "string" picture is no longer 
valid and the hole energy is minimized by a ferromagnetic 



"spin bag" , in which the hole can hop without modifying 
the background. For a circular bag with radius i? 3> a = 
1, the hole energy is given hy^ 

Eo = -Zt + ^Vli , 

where ^ = 8.53 and ^ = 5.61 for the square and the 
honeycomb lattice, respectively. We would like to point 
out, that these considerations are only valid in the case 
of the t-J model, because in the absence of quantum 
fluctuations, it is energetically favorable to align the spins 
in the plane perpendicular to the quantization axis'^^. 
The single-hole ground state therefore carries a total spin 
SL, = 0. 



4. SCBA for the t-J" model 



In the SCBA for the Ising case ©, the Bogoliubov 
parameters and the dispersion are given by Wq = 1 , Wq = 
0, and eq = J^Z/2. The Dyson equation therefore 
becomes momentum independent 



Gdie) 



1 



e - Zt^Gd{e- J^Z/2) 
and can be solved analytically'*^" 



(11) 



Gd (e) 



1 



J_ 



Ztj_ 



2e/(.7^Z)-l 



(12) 



Here J^, {x) is a Bessel function of the first kind. Note 
that the solution H12|) is not an exact solution of the 
t-J^ model. The hole, dispersionless in the SCBA, is 
mobile even in the absence of quantum fluctuations due 
to Trugman processes^. The simplest one corresponds 
to the hole motion around an elementary plaquette of 
the lattice and effectively leads to the propagation of 
the hole to a next-nearest neighbor without leaving over- 
turned spins behind. These processes are of sixth and 
tenth order in t/ for the square and the honeycomb 
lattice, respectively, and therefore negligible. More im- 
portantly, it was shown in Ref. ^3 using a diagrammatic 
study, that the inclusion of higher order terms (non- linear 
hole-magnon, static hole-magnon, direct hole-hole and 
magnon-magnon) can lead to substantially different re- 
sults. The intuitive picture emerging from that work is as 
follows: the dynamic hole-magnon interaction considered 
in the Dyson equation Hll() does not take into account the 
presence of the hole, which reduces the number of possi- 
ble sites the hole can hop to and also the energy associ- 
ated with a flipped spin. In a flrst step, the hole can hop 
to Z nearest neighbors, leaving a flipped spin with energy 
J^(Z — l)/2 behind. For the second hopping, there are 
only Z — 1 possible targets and flipping a spin costs an 
energy J^(Z — 2)/2. The Dyson equation Hll|) therefore 
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FIG. 6: Single-hole energy of the t-J' model on the square 
lattice. The SCBA results are dispersionless and based on the 
usual and modified Dyson equations ^ and respectively. 
The series expansion results are plotted for the minimum T of 
the dispersion, which is shown in the inset for J" /t = 0.3. 



has to be replaced by 

Gd (e) = 



1 



Gd (e) = 



{Z 

1 



l)/2) 



e-{Z-l)t^Gd {e-J- {Z-2)/2) 



(13) 



In this modified approach, the SCBA gives the exact re- 
sult in the limit J/t oo, where the ground state en- 
ergy is equal to = -2ZtV( J^(Z - 1)). For the fuU t-J 
model, these geometrical considerations are believed to 
be less important, because quantum fluctuations restore 
overturned spinsSS. 



IV. RESULTS FOR THE t-J^ MODEL 

Studying the simplified t-J'^ model has the advantage 
that series expansion methods are perfectly suited for this 
task and that they can therefore serve as a benchmark for 
other methods. However, it is difficult to make predic- 
tions for the t-J model on the basis of these results, since 
the modifications arising from quantum fluctuations are 
hard to quantify. Nevertheless, a comparison between 
the usual and modified SCBA on the two lat- 
tices gives interesting insight in the importance of higher 
order terms discussed in Sec. IIII C 41 and a comparison 
with series expansion results allows to estimate the in- 
fluence of Trugman processes. The common interpreta- 
tion of the hole propagation in the t-J^ model relies on 
the "string picture" : as the hole moves in the rigid Neel 
state, it scrambles the antiferromagnetic order and leaves 
a string of overturned spins behind. The cost in energy 
for these hoppings is thus proportional to the distance 
and the problem can be interpreted as a particle moving 
in a linear potential^ and the hole energies then show the 
well known (J^/t)^^^ behavior. From a numerical solu- 
tion of the modifled Dyson equation we find that 
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FIG. 7: Single-hole energy of the t-J'' model on the honeycomb 
lattice. Trugman processes are less important than on the 
square lattice, as can be seen from the remarkable agreement 
between the modfied SCBA and series expansion results, and 
also from the narrow bandwidth of the dispersion shown in the 
inset for J' /t = 0.2. 



the hole energies (of the ground state and the first few 
excited states) can be fitted by 



E/t = -2^/Z^ + a (jyt) + P {J'/t) 



2/3 



with coefficients a and /3 depending on the lattice ge- 
ometry and the state under consideration. An analytical 
estimation of /3 can be found in Ref. 



A. Square lattice 

Fig. shows the well known results for the single- hole 
ground state energies of the t-J^ model on the square 
lattice obtained via series expansionS^ and SCBA2Si2i. 
Since the SCBA does not take into account closed-loop 
paths, its results are dispersionless and we compare it 
with the series expansion results at k = F. The disper- 
sion for J^ /t = 0.3 is shown in the inset. It essentially 
represents hopping to next-nearest neighbors by loops 
turning 1.5 times around elementary plaquettes, well de- 
scribed by i?(k) = A {cos{kx + ky) -\- cos{kx — ky)) -\- C. 
It is interesting to note that in contrast to the honey- 
comb lattice, the dispersion does not include hoppings 
to all next-nearest neighbors, because only sites on the 
same plaquette can be reached within a single loop. Next- 
nearest neighbors on adjacent plaquettes require two suc- 
cessive loops and are therefore less important. As can be 
seen in Fig.|Bl the modified Dyson equation l|13|) leads to 
a very accurate description of the hole. 



B. Honeycomb lattice 

The single-hole energies of the t-J'^ model on the 
honeycomb lattice are presented in Fig. with an 
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inset showing the dispersion for /t = 0.2 ob- 
tained via series expansion. As expected, the band- 
width is much smaUer than on the square lattice, 
because Trugman processes are of higher order and 
therefore less important. Similar to the square 
lattice, the dispersion represents hoppings to next- 
nearest neighbors and can therefore be described by 
i?(k) = A{2 cos (VSfcj,) -f 4cos (\/3/2fcy) cos (3/2fc:r)) + 
C. Because all next-nearest neighbors can be reached 
by one elementary loop, the dispersion corresponds to 
hopping inside a given sublattice. Results based on the 
modified Dyson equation agree well with series ex- 
pansion calculations. However, it is interesting to note 
that the corrections due to the presence of the hole are 
more important for the honeycomb than for the square 
lattice. For instance, for /t = 0.1 the difference for the 
honeycomb lattice is 12%, but for the square lattice it is 
only 4%. 



V. RESULTS FOR THE t-J MODEL 

The t-J model on the square lattice has been the object 
of numerous studies, involving many different methods. 
In order to compare our new findings for the honeycomb 
lattice with the well known calculations for the square 
lattice, we briefly review some of these previous results 
and in addition, present a new calculation of the quasi- 
particle residue using series expansion. 



A. Square lattice 

Fig. IHta) shows the contour plot of the single-hole 
dispersion of the t-J model on the square lattice for 
J /t = 0.3. One clearly identifies four hole pockets cen- 
tred at S — (±^,±-|), which are stretched along the 
face of the magnetic Brillouin zone. The maximum is 
found at r = (0,0). The dispersion along FMST is 
shown in Fig. |Hl[b). Apart from a global offset, series 
expansion and SCBA agree perfectly well. The main 
features, namely the dispersion minimum at S and the 
flat band along the edge of the magnetic Brillouin zone 
already emerge by applying second order perturbation 
theory, and can therefore be described by hopping to 
next-nearest neighbors, i.e. hopping inside a given sub- 
lattice. Above J/t K, 0.4 the dispersion can thus be 
fitted by E (k) = 2A (cos (fc^. + ky) + cos [k^ ~ ky)) + 
2B{cos{2kx) -\-cos{2ky)) + C. Compared to the t-J^ 
model, where only hopping to next-nearest neighbors 
along the diagonals are relevant and hence B ~ 0, the dis- 
persion extrema are shifted from M to S and the overall 
shape changes substantially. For smaller values of J/t, 
the hole can hop further and additional hopping terms 
have to be taken into account. These contribution lead 
to more pronounced hole pockets. The agreement be- 
tween the SCBA and series expansion is very good up 
to J/t « 0.7, as can be seen in Fig. |Htc), showing the 



bandwidth as a function of J/t. The static magnon-hole 
interaction becomes more important with increasing J/t, 
and its neglection leads to the observed underestimation 
of the bandwidth in the SCBA. The residues Zc (k) for 
the t-J model on the square lattice for J/t — 0.3 are 
shown in Fig. |Htd). The main characteristic obtained 
via series expansion and SCBA is the almost complete 
suppression of the quasi-particle weight at X. 

B. Honeyomb lattice 

Let us start this paragraph by examining the single- 
hole spectral functions for the honeycomb lattice shown 
in Fig. |5| For the square lattice, extensive discussions 
on this topic can be found in Refs.|23 andl2ll Through- 
out the hole Brillouin zone, there is a well-defined quasi- 
particle ground state. Apart from a global energy shift, 
ED and the SCBA lead to similar results concerning the 
weight and the shape of this lowest-lying peak. How- 
ever, the structure of the excited states is different. The 
SCBA predicts several well pronounced quasi-particle 
excitations (which can be explained according to the 
"string-picture" developed for the t-J^ model), whereas 
ED results show a rather featureless incoherent tail for 
k = K, L,W and a second dominant peak at much higher 
energies at k = F, L. Given the simplicity of the SCBA, 
the overall agreement of the spectral functions is in our 
opinion quite remarkable, especially when focussing on 
ground state properties. 

Fig. llOf a) shows a contour plot of the single-hole dis- 
persion for J/t = 0.2. The circular hole pockets are lo- 
cated in the corners of the magnetic Brillouin zone W = 
(±2f ,±^) and W = (0,±^), and the maximum is 

found at F = (0, 0). The dispersion along the irreducible 
wedge TKWr is shown in Fig. llOf b). Similar to the 
square lattice, one finds a very fiat band along the edge 
of the Brillouin zone. Above J/t « 0.7, the dispersion is 
well described by next-nearest neighbor hopping i?(k) — 
A (2 cos (VSky) + 4cos {V3/2ky) cos (3/2fc^)) + C. Be- 
low J/t w 0.7 further hopping terms become increas- 
ingly important and lead to more pronounced hole pock- 
ets, analoguos to the square lattice. As can be seen in 
Fig. llOf b). the overall shape of the dispersion is very sim- 
ilar in ED calculations, the SCBA and series expansion. 
However, there is a discrepancy concerning the band- 
width: while ED (with up to 32 sites) and the SCBA 
lead to quantitatively similar results for J/t < 0.3, series 
expansion predicts a substantially larger bandwidth, see 
Fig lior c). From a general point of view, the validity of 
the SCBA is limited to small J/t, where the neglected 
hole-magnon interactions are small. The results for the 
t- model imply that these corrections are more impor- 
tant for the honeycomb lattice than for the square lattice, 
because of the reduced number of nearest neighbors. This 
implication is confirmed by the good agreement with ED 
calculations for J/t < 0.3. Still, the noticeably larger 
bandwidth obtained by series expansion is quite puz- 
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FIG. 8: Single-hole properties of the t-J model on the square lattice: (a) Contour plot of the dispersion for J/t — 0.3, with 
maximum atV = (0,0) and minima (hole pockets) at S — (±-|,±-|). (b) Dispersion for J/t = 0.3 along the irreducible wedge 
of the Brillouin zone, see Fig.\^for notations. Apart from a global offset, the SCBA and series expansion results are in good 
agreement, (c) Bandwidth W/t as a function of J/t. The agreement between series expansion and SCBA is very good up to 
J/t ~ 0.7. (d) Quasi-particle residue Zc for J/t = 0.3. 




K to W practically coincides with ED and the SCBA. 
Interestingly, the three methods lead to almost identical 
results for the quasi-particle weight shown in Fig. Iiur d) . 
Analogous to the square lattice, the weight is substan- 
tially reduced but still finite at L. In view of the notice- 
ably different predictions concerning the bandwidth, it 
would be interesting to compare these results with quan- 
tum Monte-Carlo calculations, which could elucidate the 
problem from a different point of view. 



VI. CONCLUSION 



FIG. 9: Comparison between single-hole spectral functions on 
the honeycomb lattice obtained with the SCBA and exact di- 
agonalization calculations for J/t = 0.2. 



zling. It seems that the difference is essentially due to 
the very steep band crossing the Brillouin zone from T 
to K and W , because the flat part along the edge from 



For the first time, we have calculated single-hole prop- 
erties of both the t-J^ and the t-J model on the honey- 
comb lattice using exact diagonalization, series expansion 
methods and the self-consistent Born approximation. We 
restrict our analysis to the physically relevant regime 
J/t < 1. Analogous to calculations for the very similar 
square lattice, we find an almost flat quasi-particle band 
along the edge of the magnetic Brillouin zone and well 



10 




1.2 



^ 0.6 



0.4 
0.2 



1 

<^ Series 


1 1 1 


--SCBA 




a-aED (18 sites) 




v- vED (24 sites) 








/ ^ A 




1 


^'^^ - 







0.2 



0.4 0.6 
J/t 



0.8 1 




- Series 
SCBA 
-ED (24 sites) 

I (d) 



w 



FIG. 10: Single-hole properties of the t-J model on the honeycomb lattice: (a) Contour plot of dispersion for J/t = 0.2. 
The maximum is found at T = (0, 0) and six circular hole pockets are located in the comers of the magnetic Brillouin zone 



(b) Dispersion for J/t = 0.2 along the irreducible wedge of the Brillouin zone, 



:^,±^) andW = {0,±^^, 
Fig. d for notations. The overall shape of the dispersion is the same for all three methods, and essentially due to hopping 
to next-nearest neighbors, but series expansion predicts a substantially larger bandwidth, (c) Bandwidth W/t as a function of 
J/t. The SCBA is valid only at small J/t, where it agrees well with ED. (d) Quasi-particle residue for J/t = 0.2. 



defined circular hole pockets in the corners. The spec- 
tral functions show well-defined quasi-particle peaks at 
the bottom of the spectrum throughout the whole Bril- 
louin zone. We find a good agreement between the three 
methods applied concerning the shape of the single-hole 
dispersion and the quasi-particle weight. An analysis of 
the bandwidth clearly shows that the self-consistent Born 
approximation is quite accurate for J/t < 0.3, but leads 
to a substantial underestimation of the bandwidth for 
larger J/t. Series expansion on the other hand seems 
to overestimate the bandwidth in the small J/t region. 
It is likely that this bandwidth disagreement is due to 
the reduced "dimensionality" of the honeycomb lattice, 



because the number of nearest neighbors {Z — 3) lies 
between the square lattice {Z — 4) and the linear chain 
{Z = 2). This makes the problem especially interesting 
and a Monte-Carlo calculation would be very helpful to 
finalize the findings of the present work. 
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